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1 Introduction. 



Most microscopic many-body quantum systems are subject to strong in- 
teractions which act between their constituents. In general, there exist no 
analytical methods to treat exactly strongly interacting systems apart from 
the assumption of trial wave functions able to diagonalize the Hamiltonian of 
integrable models, like the Bethe Ansatz (BA) for specific one dimensional 
systems (1-D) [1], the BCS hypothesis which explains the supraconductivity 
and others like the Neel state, the Resonant Valence Bond {RVB) spin liquid 
states proposed by Anderson |2], the Valence Bond Crystal states (VBC). . . . 
They are aimed to describe 2D systems but there remains the problem of 
their degree of realism, i.e. their ability to include the essentials of the inter- 
action in strongly interacting systems p]. The complexity of the structure 
of such systems leads to the diagonalization of the Hamiltonian numerically 
which must in general be performed in very large Hilbert space although 
the information of interest is restricted to the knowledge of a few low en- 
ergy states generally characterized by collective properties. Consequently it 
is necessary to manipulate very large matrices in order to extract a reduced 
quantity of informations. 

Non-perturbative techniques are needed. During the last decades a con- 
siderable amount of procedures relying on the renormalization group concept 
introduced by Wilson [3] have been proposed and tested. Some of them are 
specifically devised for quantum spin systems, like the Real Space Renormal- 
ization Group (RSRG) [H [6] and the Density Matrix Renormalization Group 
(DMRG) [BIB]. 

We propose here a non-perturbative approach which tackles this question 
[9]. The procedure consists of an algorithm which implements a step by step 
reduction of the size of Hilbert space by means of a projection technique. It 
relies on the renormalization concept following in spirit former work based 
on this concept [IHl [HI [12] • Since the reduction procedure does not act in 
ordinary or momentum space but in Hilbert space, it is universal in the sense 
that it works for any kind of many-body quantum system. 

The properties of physical systems can be investigated in different repre- 
sentations. In the present work which deals with frustrated spin ladders the 
common S'?7(2)-representation is confronted with the S'0(4)-representation 
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[ini [HI 13 in order to work out the energies of the low-lying levels of the 
spectrum of these systems. The efficiency of one or the other representation 
in terms of the number of relevant basis states characterizing the ground 
state wave function is tested in connection with the reduction process. 

The outline of the paper is the following. In section 2 we present the 
formal developments leading to the secular equation in the reduced Hilbert 
space. Section 3 is devoted to the application of the algorithm to frustrated 
quantum spin ladders with two legs and one spin per site. We analyze the 
outcome of the applied algorithm on systems characterized by different cou- 
pling strengths by means of numerical examples, in bases of states developed 
in the SU{2) and S'0(4)-representations and compare the results obtained 
in both cases. Conclusions and further possible investigations and develop- 
ments are presented in section 4. 



2 The reduction algorithm. 

2.1 Reduction procedure and renormalization of the 
coupling strengths. 

We consider a system described by a Hamiltonian depending on a unique 
coupling strength g which can be written as a sum of two terms 

H = Ho + gH, (1) 

The Hilbert space H^^-* of dimension is spanned by an orthonormalized 
arbitrary set of basis states i = 1, ■ ■ ■ , N}. In this basis an eigenvector 

|\&;'^^) takes the form 

N 

i^r)=E«r(^^"^)i^^) (2) 
1=1 

where the amplitudes {0'if\g^^^)} depend on the value (7^^^ of g in n^^^ 

Using the Feshbach formalism [16] the Hilbert space may be decomposed 
into subspaces by means of the projection operators P and Q, 
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In practice the subspace PTi^^^ is chosen to be of dimension dimP'H^^'' = 

— 1 by ehmination of one basis state. The projected eigenvector P\'^[^^) 
obeys the Schroedinger equation 

H.fA^r)p\ir} = Ar'pifl-r) . (4) 

where Heff{X[^^) is the effective Hamihonian which operates in the sub- 
space Pl-i'^^\ It depends on the eigenvalue A;^"* which is the eigenenergy 
corresponding to 1^1/;^'') in the initial space "H^^^ The coupling strengths 
^(^) which characterizes the Hamiltonian H'^^'^ in "H*^^^ is now aimed to be 
changed into f?*-^"^^ in such a way that the eigenvalue in the new space "H*^^"^-* 
is the same as the one in the complete space 

Aj"-' = Aj"' (5) 

The determination of g^^~^^ by means of the constraint expressed by Eq. 
is the central point of the procedure. It is the result of a renormalization 
procedure induced by the reduction of the vector space of dimension to 

— 1 which preserves the physical eigenvalue X[^^ . 

In the sequel P\'^[^^) is chosen to be projection of the ground state eigen- 
vector I^E'^^'') {I = 1) and A^^'' = A^^"^'' = Ai the corresponding eigenenergy. 
In ref. [S] it is shown how g^'^~^^ can be obtained as a solution of an alge- 
braic equation of the second degree. One gets explicitly a discrete quadratic 
equation 

where 

6(^-1) = aSf i7^^(Ar - «i) + Fi^(Ar - (8) 
c(^-^) = -aSr^(AS^^-«i)(AS^)-a^) (9) 

with 

Af-l 

i^iiv = $^aSf (<|.i|iJi|$,) 

i=l 
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N-1 
i=l 

H,, = ($.|ifi|<f,) 

and 

ai = {<l>i\Ho\<l>i), i = l,--- ,N 

The reduction procedure is then iterated in a step by step decrease of the 
dimensions of the vector space, A^i— i-A^ — 1 N — 2 ... leading at each 
step to a couphng strength g^^~''^ which can be given as the solution of 
a flow equation in a continuum limit description of the Hilbert space. The 
procedure can be generalized to Hamiltonians depending on several coupling 
constants which experience a renormalization during the reduction procedure 
under further constraints |18j . 



2.2 Outline of the reduction algorithm. 

We sketch here the different steps of the procedure. 

1— Consider a quantum system described by an Hamiltonian H^'^^ which 
acts in an A^-dimensional Hilbert space. 

2— Compute the matrix elements of the Hamiltonian matrix H^'^'^ in a 
definite basis of states = 1, . . . ,N}. The diagonal matrix elements 
{ej = are arranged either in increasing order with respect to 
the {ei} or in decreasing order of the absolute values of the ground state wave 
function amplitudes \a'ii \g^^'^)\ [T7] . 

3— Use the Lanczos technique to determine A^^-* and \'^[^\g^^'^)) [T91I20] . 

4— Fix g^^~^^ as described in section 2.1. Take the solution of the alge- 
braic second order equation closest to c/(^) (see Eq.®). 
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5— Construct H^^ = Hq + g^^ ^^Hi by elimination of the matrix ele- 
ments of H'^^'^ involving the state \^n)- 

6— Repeat procedures 3, 4 and 5 by fixing at each step A; {A; = 1, . . . , A^ — 
Nmin], = A^^^ = Ai- The iterations are stopped at Nmin correspond- 
ing to the limit of space dimensions for which the spectrum gets unstable. 

2.3 Some remarks. 

• The procedure is aimed to generate the energies of the low-energy ex- 
cited states of strongly interacting systems and possibly the calculation 
of further physical quantities. 

• The implementation of the reduction procedure asks for the knowl- 
edge of Ai and the corresponding eigenvector '^■*) at each step k 
of the reduction process. The eigenvalue Ai is chosen as the physical 
ground state energy of the system. Eigenvalue and eigenvector can 
be obtained by means of the Lanczos algorithm [SI HSl 120] which is 
particularly well adapted to very large vector space dimensions. The 
algorithm fixes A^^ = A^^"* and determines at each step. 

• The process does not guarantee a rigorous stability of the eigenvalue 
Ai. I^I/i^ ^ which is the eigenvector in the space and 
the projected state P|^[^"*-'^) of into n^^-^~^) may differ 
from each other. As a consequence it may not be possible to keep 
y{N k 1) j^jgQj^Q^giy equal to A^^ ^'^ = Ai. In practice the degree of 
accuracy depends on the relative size of the eliminated amplitudes 
{^i(N~k)i9^'^~^^)}- '^^^^ point will be tested by means of numerical 
estimations and further discussed below. 

• The Hamiltonians of the considered ladder systems are characterized 
by a fixed total magnetic magnetization Mfot- We work in subspaces 
which correspond to fixed Mfot- The total spin Stot is also a good 
quantum number which defines smaller subspaces for fixed Mtot pT] . 
We do not introduce them here because projection procedures on Stot 
are time consuming. Furthermore we want to test the algorithm in 
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Figure 1: Top: the original spin ladder. The coupling strengths are indicated 
as given in the text. Bottom: The ladder in the SO(4)-representation. See 
the text. 

large enough spaces although not necessarily the largest possible ones 
in this preliminary tests considered here. 

3 Application to frustrated two-leg quantum 
spin ladders. 

3.1 The model 

3.1.1 SU(2)-representation. 

Consider spin- 1/2 ladders ^21123] described by Hamiltonians of the following 
type and shown in Fig.l 

L 

i=l <ij> <ij> (ij) 

+ J2C Si^Sj^ 

The indices 1 or 2 label the spin 1/2 vector operators Si^ acting on the 
sites i on both ends of a rung, in the second and third term i and j label 
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nearest neighbours, here j = i + 1 along the legs of the ladder. The fourth 
and fifth term correspond to diagonal interactions between sites located on 
different legs, j = i + 1. L is the number of sites on a leg (Fig. 1) where 
Jic = J2c = Jc- The coupling strengths Jj, J;, Jc are positive. 

As stated above the renormalization is restricted to a unique coupling 
strength, see Eq. ([T]). It is implemented here by putting ifg = and H^^^ = 
g^Hi where g^^^ = Jt and 

L 

4=1 <ij> <ij> 

where jti = Ji/Jt, 7c = Jc/ Jt- These quantities are kept constant and 
gW = y^iii \)Q subject to renormalization during the reduction process. 

One should point out that the renormalization does not change if one 
chooses another coupling parameter as a renormalizable parameter, here Ji 
or Jc, because they are related to each other at the beginning of the reduction 
procedure by the ratios jti = Ji/Jt and 7c = Jc/ Jt- 

The basis of states {|$p) ,p = 1, • • • , N} is chosen as 

2L 

|$p) = |l/2 mi,. ..,1/2 mi, ...,1/2 maL, m^ = M^^ = 0) 

i=l 

with {rrii = +1/2,-1/2}. 

3.1.2 SO(4)-representation. 

Different choices of bases may induce a more or less efficient reduction pro- 
cedure depending on the strength of the coupling constants Jt, Ji, Jc- This 
point is investigated here by choosing also a basis of states which is written 
in an 5'0(4)-representation. 

We replace (sji,Si2) corresponding to dimers by {Si,Ri). By means of a 
spin rotation [HI [15] 

Si, = ^{S, + Ri) . (12) 
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The Hamiltonian, Eq. f llOp . can be expressed in the form 



(13) 



4=1 <ij> <ij> 

The structure of the corresponding system is shown in the lower part of Fig. 
1. Here Ji = (J; + Jc)/2, J2 = (Ji — Jc)/'^ and as before Jic = J2C = Jc- 
The components S^i^\ S\^\ S\'''' and R^^\ B!f\ Blf"^ of the vector operators 
Si and Ri are the SO {4) group generators and < ij > denotes nearest neigh- 
bour indices. 

In this representation the states {\SiMi)} are defined as 

\SiMi) = J2 (1/2 mi 1/2 m2\SiM,)\l/2 mi)i|l/2 ms)^ 

mi ,m2 

along a rung are coupled to S'j = or S'j = 1. Spectra are constructed in this 
representation as well as in the S'?7(2)-representation and the states {|$p)} 
take the form 

L 

1$^) = l^iMi, S,Mi, SlMl, 5Z M, = Mtot = 0) 

i=l 

3.2 Test observables 

In order to quantify the accuracy of the procedure we introduce different 
test quantities in order to estimate deviations between ground state and low 
excited state energies in Hilbert spaces of different dimensions. The stability 
of low-lying states can be estimated by means of 

/ (N) _ (ri)x 

p(i) = \^ — X 100 with i = l,...,4 (15) 

where e^"^ = /2L with n = {N — k) corresponds to the energy per site at 
the ith physical state starting from the ground state at the kth iteration in 
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Hilbert space. This quantity provides a percentage of loss of accuracy of the 
eigenenergies in the different reduced spaces. 

A global characterization of the ground state wavefunction in different 
representations can also be given by the entropy per site in a space of dimen- 
sion n 



s = -^5^P,/nP, with i^, = |($(")|M/("))|2 = |a;:)|2 (16) 

i=l 

which works as a global measure of the distribution of the amplitudes 
in the physical ground state. 

In the remaining part we work out the spectra of different systems and 
compare results obtained in the two representations introduced above. 

3.3 Spectra in the SU(2)-representation. 

Results obtained with an S'f/(2)-representation basis of states are shown in 
Figs.{2 - 3). 

3.3.1 First case: L= 6, Jt=15, J/=5, Jc=3 

We choose the basis states in the framework of the M-scheme corresponding 
to subspaces with fixed values of the total projection of the spin of the {|$j)}, 
Mtot = 0. 

In the present case Jt > Ji, Jc- The dimension of the subspace is reduced 
step by step as explained above starting from N = 924. As stated in section 
2.2 the basis states are ordered with increasing energy of their diago- 

nal matrix elements {ej} and eliminated starting from the state with largest 
energy e^v. 

Deviations of the energies of the ground and first excited states from their 
initial values at = 924 can be seen in Figs.([2^-b) where the p{i)^s defined 
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above represent these deviations in terms of percentages. As seen in Fig. 
the ground state of the system stays stable down to n ~ 50 where n is the 
dimension of the reduced space. The couphng constant Jt does not move 
either down to n ~ 300, see Fig.(|2l:;). Figs.([2^-b) show the evolution of the 
first excited states which follows the same trend as the ground state. 

For n < 50 the spectrum gets unstable, the renormalization of the cou- 
pling constant can no longer correct for the energy of the lowest state. Indeed 
the coupling constant Jt increases drastically as seen in Fig.(|2l:;). The reason 
for this behaviour can be found in the fact that at this stage the algorithm 
eliminates states which have an essential component in the state of lowest 
energy. The same message can be read on Fig.([21i), the drop in the entropy 
per site s is due to the elimination of sizable amplitudes {aii}. 

3.3.2 Second case: L= 6, Jt=5.5, Ji=5, Jc=3 

Contrary to the former case the coupling constant Jt along rungs is now of the 
order of strength as Ji, Jc- Results are shown in Fig.(|3]). The lowest energy 
state is now stable down to n ~ 100. This is also reflected in the behaviour 
of the excited states which move appreciably for n < 200. Fig.(|3t) shows 
that the coupling constant Jt starts to increase sharply between n = 300 and 
n = 200. It is able to stabilize the excited states down to about n = 200 and 
the ground state down to n = 70. The instability for n < 70 reflects in the 
evolution of the p(i)'s, Figs.(|3^-b) which get of the order of a few percent. 
The entropy Fig.(|2Ji) follows the same trend. 

Comparing the two cases above and particularly the entropies Fig. ([311) 
and Fig.([2ji) one sees that the stronger Jt the more the amplitude strength 
of the ground state wavefunction is concentrated in a smaller number of basis 
state components. The elimination of sizable components of the wavefunc- 
tion leads to deviations which can be controlled down to a certain limit by 
means of the renormalization of Jt. One sees that large values of Jt favour a 
low number of significative components in the low energy part of the spec- 
trum in a SU (2)- representation. 
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3.3.3 Remark 



In Figs.(|2^-b) it is seen that the ground and first excited states show "bunches" 
of energy fluctuations. In the ground state the peaks are intermittent, they 
appear and disappear during the space dimension reduction process. They 
are small in the case where Jt = 15 but can grow with decreasing Jf as it 
can be observed in Figs.(|3^-b) for Jt = 5.5. The subsequent stabilization of 
the ground state energy following such a bunch shows the effectiveness of the 
coupling constant renormalization which acts in a progressively reduced and 
hence incomplete basis of states. 

These bunches of fluctuations are correlated with the change of the num- 
ber of relevant amplitudes (i.e. amplitudes larger than some value e as ex- 
plained in the caption of Fig.(|6])) during the reduction process. 

Consider first the case where Jt > Ji, Jc- One notices in the caption of 
Fig. (jn^) that down to n ~ 300 the number of relevant amplitudes defined in 
Fig.dS]) stays stable like the ratios {p{i)} in Figs.(I2^-b). For 158 < n < 300 
these ratios change quickly. A bunch of fluctuations appears in this domain of 
values of n as seen in Figs. ([2^-b) and correspondingly the number of relevant 
amplitudes decreases steeply. For 60 < n < 158 the ratios {p{i)} stay again 
stable as well as the number of relevant amplitudes. The {p{i)} in Figs.([2^-b) 
almost decrease back to their initial values. The analysis shows that these 
bunches of fluctuations signal the local elimination of relevant contributions 
of basis states to the physical states in the spectrum. The stabilization of 
the spectra which follows during the elimination process shows that renor- 
malization is able to cure these effects. 

In the case where Jt ~ J/, Jc shown in Fig.(|6]D) the relevant and irrelevant 
amplitudes move continuously during the reduction process and the corre- 
sponding {p{i)} do no longer decrease to the values they showed before the 
appearance of a bunch of energy fluctuations as seen in Figs.(|3^-b). It signals 
the fact that the coupling renormalization is no longer able to compensate 
for the reduction of the Hilbert space dimensions. 

One should mention that the evolution of the spectrum depends on the 
initial size of Hilbert space. The larger the initial space the larger the ratio 
between the initial number of states and the number of states corresponding 
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to the limit of stability of the spectrum, see ref. [IB] for explicit numerical 
examples. 

3.4 Spectra in the SO(4)-representation 

The reduction algorithm is now applied to the system described by the Hamil- 
tonian given by Eq. f lT^ with a basis of states written in the 50(4)- 

representation. Like above we consider two cases corresponding to large and 
close values of Jt relative to the strengths of the other coupling parameters. 

3.4.1 Reduction test for L = 6, Jt = 15, 5.5 and Ji=5, Jc=3 

Figs. (jUlS]) show the behaviour of the spectrum for a system of size L = 6. 
A large value of Jt, {Jt = 15), favours the dimer structure along rungs in 
the lowest energy state and stabilizes the spectrum down to small Hilbert 
space dimensions. This effect is clearly seen in Fig.(|l^), the ground state is 
very stable. The excited states are more affected, see Fig. (jib), although they 
do not move significantly. The renormalization of the coupling strength Jt 
starts to work for n ^ 50, Fig. (jit). 

The situation changes progressively with decreasing values of Jt. Figs. ([5]) 
show the case where Jt = 5.5. The ground state energy experiences sizable 
bunches of fluctuations like in the SU (2)-representation, but much stronger 
than in this last case. The same is true for the excited states which is reflected 
through all the quantities shown in Figs.([S^-b), in particular Jt, Fig.dSt). The 
entropy Figs. (jUi jSji) follows the same trend like in the 5'?7(2)-representation 
by changing the values of Jt from 15 to 5.5. The analysis of the bunches of 
energy fluctuations through the number of relevant-irrelevant amplitudes is 
shown for Jt = 15, J/ = 5, Jc = 3 in Fig. (^.l-c.2) and for Jt = 5.5, J/ = 5, 
Jc = 3 in Fig. (p). 

The results show that the renormalization procedure is quite sensitive 
to the representation chosen in Hilbert space. It is expected that essential 
components of the ground state wavefunction get eliminated early during 
the process when the rung coupling gets of the order of magnitude or smaller 
than the other coupling strengths. 
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By comparing the two representations in Fig([S]), one notices that the 
stabihty of the low-energy properties of the system in the reduced Hilbert 
space is well characterized by the number of relevant-irrelevant amplitudes 
as defined in the caption of Fig.(|6]) and the distribution of the amplitudes in 
Hilbert space, see ref. [T7]. 



3.5 Summary 

The present results lead to two correlated remarks. The efficiency of the algo- 
rithm is different in different sectors of the coupling parameter space. In the 
case of the frustrated ladders considered here the algorithm is the more effi- 
cient the stronger the coupling between rung sites Jj. Second, this behaviour 
is strongly related to the representation in which the basis of states is defined. 
The S'f/(2)-representation leads to a structure of the wave functions (i.e. the 
size of the amplitudes of the basis states) which is very different from the one 
obtained in the S'0(4)-representation. For large values of Jt the spectrum is 
more stable in the 5'0(4)-representation. For small values of Jt the stability 
is better realized in the SU (2)- representation. Finally, in the regime where 
Jt > Ji, Jc, one observes that the reduction procedure is the more efficient 
the closer J; to Jc- This effect can be understood and related to previous 
analytical work in the S'0(4)-representation [25] . 



4 Conclusions and outlook. 

In the present work we tested and analysed the outcome of an algorithm 
which aims to reduce the dimensions of the Hilbert space of states describing 
strongly interacting systems. The reduction is compensated by the renor- 
malization of the coupling strengths which enter the Hamiltonians of the 
systems. By construction the algorithm works in any space dimension and 
may be applied to the study of any microscopic A^-body quantum system. 
The robustness of the algorithm has been tested on frustrated quantum spin 
ladders. 

The analysis of the numerical results obtained in applications to quantum 
spin ladders leads to the following conclusions. 
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The present numerical applications are essentially realistic tests of 
the renormalization algorithm restricted to rather small Hilbert spaces 
which do not necessarily necessitate the use of a Lanczos procedure. We 
introduced this procedure in order to be able to extend our present work 
to much larger systems for which ordinary full diagonalization cannot 
be performed. We tested it sucessfuUy by comparing the outcome with 
that of complete diagonalization. Preliminary extended calculations 
have also been performed on higher dimensional ladders. We expect to 
present and discuss the physics of their application to frustrated sys- 
tems in further work. 

The stability of the low-lying states of the spectrum in the course of 
the reduction procedure depends on the relative values of the coupling 
strengths. The ladder favours a dimer structure along the rungs, i.e. 
stability is the better the larger the transverse coupling strength J^. 
This is the reason why the 5'0(4)-rcprcscntation is favoured when com- 
pared to the SU{2)- representation in this case. It leads to a more 
efficient basis of states as commented below (see next paragraph). 

The efficiency of the reduction procedure depends on the representation 
frame in which the basis of states is defined. It appears clearly that the 
evolution of the spectrum described in an SU{2)- representation is sig- 
nificantly different from the evolution in an S'0(4)-representation. This 
is understandable since different representations partition Hilbert space 
in different ways and favour one or the other representation depend- 
ing on the relative strengths of the coupling constants. It is always 
appropriate to work in a basis of states whose symmetry properties 
are closest to the symmetry properties of the physical system so that 
the contributions of the non-diagonal couplings are the smallest. The 
choice of an optimal basis may not always be evident. In the case of 
the ladder we consider here the 5'0(4)-representation. 

Local spectral instabilities appearing in the course of the reduction pro- 
cedure are correlated with the elimination of basis states with sizable 
amplitudes in the ground state wavefunction. One or another repre- 
sentation can be more efficient on the reduction process for a given 
set of coupling parameters because it leads to physical states in which 
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the weight on the basis states is concentrated in a different number of 
components. This point is strongly related to the correlation between 
quantum entanglement and symmetry properties which have been un- 
der intensive scrutiny, see f.i. [26] and refs. therein. 

Further points are worthwhile to be investigated: 

• Extension to larger ladders and systems of higher space dimensions 
with the help of more sophisticated numerical algorithms ES]- 

• Extension of the renormalization procedure to systems at finite temper- 
ature and more than one coupling constant renormalization [18j. 
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Figure 2: SU (2) — representation, n is Hilbert space dimension. L = 6 sites 
along a leg. Jt = 15, Ji = 5, Jc = 3 
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Figure 3: SU (2) — representation, n is Hilbert space dimension. L = 6 sites 
along a leg. Jt = 5.5, J/ = 5, Jc = 3 
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Figure 4: SO {4) — representation, n is Hilbert space dimension. L = 6 sites 
along the chain. Jt = 15, Ji = 5, Jc = 3 
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Figure 5: 50(4) — representation, n is Hilbert space dimension. L = 6 sites 
along the chain. Jt = 5.5, Ji = 5, Jc = 3 
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Figure 6: n is the dimension of the Hilbert space. Amplitudes show the 
number of relevant -irrelevant amplitudes in the ground state eigenfunction 
(relevant-irrelevant amplit in the figure). Relevant amplitudes are those for 
which {|aij| > e, (here e = 10~^), i = 1, . . . ,n}. The number of sites along 
a leg is L = 6. SU{2) — representation: (a) corresponds to Jt = 15, J/ = 5, 
Jc = 3. (b) corresponds to Jt = 5.5, Ji = 5, Jc = 3. 5*0(4) — representation: 
(c.l) and (c.2) correspond to Jt = 15, Ji = 5, Jc = 3. (d) corresponds to 
Jt = 5.5, Ji = 5, Jc = 3. 
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